Sodium-terminated zigzag graphene nanoribbon: A one-dimensional semimetal with a 

tilted Dirac cone 
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Using first principles pseudopotential density functional theory calculations, we find that terminat- 
ing zigzag graphene nanoribbons (ZGNR) with monovalent alkali atoms at a reduced concentration 
has a dramatic impact on their properties. In particular, using sodium atoms for the saturation 
of ZGNR edges at half the concentration of edge-carbon atoms make it a one dimensional, perfect 
semimetal, where the valance and conduction bands meet at only a single, Dirac-like point. Unlike 
pristine graphene, the Dirac- "cones" of Na-ZGNR is not symmetric with respect to wave vector, 
but rather it is tilted. 
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After its realization as a stable monolayer structure [l| , 
graphene has become a material of great interest for re- 
search in material science due to its impressivephysical 
properties and the potential to tune them [246|- Pris- 
tine graphene is a zero-bandgap (perfect) semimetal. The 
conduction and valance bands of graphene meet linearly 
at six distinct points (Dirac-points) of first Brillouin zone, 
and are in the form of symmetric cones at the Fermi level 
01 • Long and thin strips of graphene, known as graphene 
nanoribbons (GNR), are semiconducting and they have 
been experimentally realized with smooth edges and a 
range of widths down to sub 10 nm scale As a 

special class, zigzag-edged GNRs (ZGNR) have width 
dependent energy band gaps, and are also expected to 
present edge localized spin states [IMl. These local- 
ized spin states are ferromagnetically ordered at each 
edge and they are antiferromagnetically coupled to each 
other m. 

In computational studies of ZGNRs, to facilitate me- 
chanical stability of the structure, the edges are usually 
terminated with hydrogen atoms, which lowers the val- 
ues of its energy band gap and the induced magnetic 
moments of the edges. To manipulate the electronic 
properties of ZGNRs, modification 16|, [l7| and func- 
tionalization [18[ of edges and producing defects, dop- 
ing [H-ii] and absorbtion 1231-1251 of atoms or molecules 



on ZGNRs have been extensively studied. For instance, 
recent studies 2^ 27 1 have shown ways to obtain half- 
metaUic ZGNRs. 

In this letter we report that edge saturation of ZGNRs 
with Na atoms at a specific concentration makes them a 
zero-gap one dimensional semimetal with an asymmetric 
linear energy dispersion around the Fermi level. That 
corresponds to a Dirac cone similar to pristine two di- 
mensional (2D) graphene 0], however, with a tilted form. 
Similar to hydrogenization of the edges, Na atoms satu- 
rates the dangling bonds of edge carbon atoms, and mod- 
ify the local magnetic moment magnitudes and change 
the electronic band structure. 

Geometry optimizations of bare and Na-terminated 




FIG. 1: (Color online) Schematic representation of various 
binding geometries of Na atoms to a ZGNR edge. From left 
to right: top in-plane (TI), bridge off-plane (BO), bridge in- 
plane (BI). 



ZGNRs and calculations of their magnetic and electric 
properties were performed using the SIESTA package [28 1 
based on density functional theory (DFT) '29\. We have 
used generalized gradient approximation (GGA) for the 
exchange and correlation potential as parameterized by 
Perdew, Burke and Ernzerhof [s^. For geometry opti- 
mizations, a local relaxation has been performed using 
the conjugate gradient algorithm and the convergence 
criteria of 0.04 eVxA~^ and 10~^ eV for the forces and 
total energies, respectively, in the self-consistency cycles. 
The electrostatic potentials were determined on a real- 
space grid with a mesh cutoff energy of 300 Ry. We 
make use of non-conserving TrouUier Martins pseudopo- 
tentials (sij in the Kleinman Bylander factorized form 
32| and a double- C polarized basis set composed of nu- 
merical atomic orbitals of finite range. The Brillouin 
zone has been sampled with (1,1,70) points within the 
Monkhorst-Pack k-point sampling scheme. 

We have considered edge termination of n-ZGNRs (n 
denoting the number of zigzag rows in the ribbon) with 
monovalent atoms, i.e. H, Li, Na and K. Here we will re- 
port our results using Na terminated 9- ZGNR as a rep- 
resentative case. In FiglU different binding geometries 
of Na atoms to ZGNR determined after geometry opti- 
mizations are shown; these are, on top of an edge C atom, 
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FIG. 2: (Color online) Relaxed geometric structure and iso- 
surfaces of charge density difference of spin-up (f) and spin- 
down (D states for 9-ZGNR; bare (a), and Na terminated at 
HC (b) and FC (c) concentrations. For all geometries, pos- 
itive and negative values of the charge density difference is 
shown by red (dark) and blue (light) regions, respectively, for 
the same isosurface value of ±0.0025 electrons/A^. 



in-planc vifith the ribbon (TI); bridge site of two edge C 
atoms, cither off-plane (BO), or in plane (BI). Unlike H- 
termination, Na atoms prefer to bind at the bridge sites 
rather than the TI site; binding energy of a single Na 
atom at a BI site is 0.46 eV larger than that at a TI site. 
We considered two different concentrations of Na for sat- 
uration of the edges, full coverage (FC) where there are 
as many Na atoms as the edge C atoms, and half cover- 
age (HC) where the Na atoms are bound to every other 
bridge site only (see Figl2]for the unit cell of the ribbon 
and Na binding sites). In the HC case, equilibrium posi- 
tion of Na atoms is BI with a Na-C bond length of 2.33 
A. In the FC case, however, due to their larger atomic 
size than H, neighboring Na atoms relax towards oppo- 
site directions at the BO sites. In this case, Na-C bond 
length and vertical distance of Na to the plane of ribbon 
are 2.35 A and 1.23 A respectively. Accordingly, the 
binding energy per Na atom reduces to 2.84 eV at FC 
from 2.91 eV at HC. 

We calculated spin-dependent charge densities of the 
optimized structures and their difference (Ap — p^ — p^). 
The characteristic edge-localized spin states of ZGNRs 
are preserved for both HC and FC concentrations of Na 
terminated 9-ZGNR (see Figl2]). Using MuUiken popula- 
tion analysis we have determined spin dependent atomic 
charges, total valance charge and the spin magnetic mo- 
ments of the atoms. Total magnetic moments of the 
ribbons are found zero for all cases. As seen in the 
figure, there are no net induced magnetic moments on 
the Na atoms either. The penetration of the spin imbal- 
ance on the C atoms, Ap, toward the center of ribbons 
is considerably shorter for the HC case. Binding of Na 
atoms reduces the magnetic moments of the edge carbon 



atoms from fi = 1.21 /is to /i = Q.llfxs for HC, and 
to /i = 0.28 /is for FC. This magnetic moment reduc- 
tion is comparable to the full hydrogenation of ZGNRs 
where edge carbon magnetic moments were reported as 



^l = 0.26 /^B [331. When a Na atom binds to ZGNR 



edge, its valance charge is partially transferred to the 
ZGNR due to lower electronegativity of Na than C. At 
HC (FC) concentration, Na transfers 10% (15%) of its to- 
tal valance charge. Correspondingly, total valance charge 
of each edge carbon atom of 9-ZGNR increases from 4.06 
for bare ZGNR to 4.13 at HC and 4.20 at FC concentra- 
tion. 

In Fig. [3] we display the energy band diagram and den- 
sity of states of Na terminated 9-ZGNR at HC and FC 
concentrations. Each case modifies the electronic prop- 
erties of the ribbon differently. We calculate that bare 
9-ZGNR is a direct gap semiconductor with Eg = 0.70 
eV. Fully hydrogenated 9-ZGNR has a lower band gap 



of i?g = 0.50 eV, consistent with a previous report [34 



Similarly, 9-ZGNR terminated with Na at FC concen- 
tration has a reduced and slightly indirect band gap of 
Eg = 0.44 eV. 

At HC concentration of Na termination the 9-ZGNR 
becomes a zero-gap semimetal with an interesting band 
structure. The valance and conduction bands intersect 
at a single point. In Fig.lSja), the Dirac-like point is 
located at fco — 0.469 A"-'^ in the first Brillouin zone 
{kz = 0.625 A-i). The tilted Dirac cone is formed by 
a steep linear band crossing a relatively flat band at the 
Fermi level. These bands are mainly derived from car- 
bon px and combination of carbon py and Na s orbitals, 
respectively. Since the present ribbon structure is one di- 
mensional, the density of states (DOS) at the Fermi level 
does not vanish unlike 2D-graphene, and also, singulari- 
ties characteristic to one-dimensional systems is evident 
in the DOS plot. 

Energy band structures showing Dirac cones are quite 
rare. 2D-graphene is known for having symmetric and 
isotropic Dirac cones at corners of the first Brillouin zone. 
In this context, recently another two dimensional struc- 
ture, which is an organic compound, has been reported 
as a unique material having a gapless band diagram and 
asymmetrical linear energy dispersion 



35|. Na termi- 



nated ZGNR at HC concentration that we present in this 
study is a one dimensional member of this class of mate- 
rials having tilted Dirac "cones" . 

Around the Fermi level, E vs k relation can be approx- 
imated in the form. 



Ex{k) = wq {k - ko) + Xw\k — ko\ 



(1) 



where A plays the role of band index (+1 for conduc- 
tion, and —1 for valance band), w — 2.517 eVxA and 
wq — —2.277 eVxA are effective velocities. In this rep- 
resentation the Wo parameter determines the degree of 
tilting of the Dirac cone. The slopes of the bands at the 
Fermi level are wq+w = 0.240 eVx A and wq — w = 4.794 
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FIG. 3: (Color online) Band structures and density of states 
diagrams for 9-ZGNR. (a) HC, and (b) FC concentration of 
Na. The inset magnifies the tilted Dirac cone at the Fermi 
level for the HC case. 



eVxA, so that an asymmetry of about a factor of 20 
in the slopes of the fermion and antifermion bands is 
present. 

We find that the "tilted Dirac-cone" shaped band 
structure is quite robust with respect to the ribbon width 
and the approximations involved, such as exchange cor- 
relation functional being LDA or GGA. Using Li or K 
instead of Na also produces similar electronic structures, 
however, the valance and conduction bands intersect at 
slightly below or above the Fermi energies, respectively. 

In summary, we have performed first-principles to- 
tal energy calculations to investigate edge saturation of 
ZGNRs with Na atoms. Binding a Na atom to every 
other bridge site of the zigzag edge leads to formation 
of linear bands crossing at the Fermi level. To our 
knowledge, such tilted Dirac "cones" in a one dimen- 
sional material has not been reported before [s^, and 
may find interesting spintronics applications. Such a one- 
dimensional semi-metal with charge carriers having high 
mobility is desirable for realization of nanoscale devices. 
Our preliminary results also show that edge saturation of 
ZGNRs with other monovalent alkali metals (Li and K) 
have similar effects, but those will be reported elsewhere. 

This work was partially supported by IzTech-BAP pro- 
gram through project no: 2010IYTE26. 
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